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<N i Abstract 

We investigate the nature of randomness in disordered packings of frictional spheres. We cal- 
culate the entropy of 3D packings through the force and volume ensemble of jammed matter, a 
mesoscopic ensemble and numerical simulations using volume fluctuation analysis and graph theo- 
retical methods. Equations of state are obtained relating entropy, volume fraction and compactivity 
§ \ characterizing the different states of jammed matter. At the mesoscopic level the entropy vanishes 

at random close packing, while the microscopic states contribute to a finite entropy. The entropy 

£ 

of the jammed system reveals that the random loose packings are more disordered than random 
close packings, allowing for an unambiguous interpretation of both limits. 
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Filling containers with balls is one of the oldest mathematical physical puzzles known to 
scientists [l]. Apart from its mathematical significance, this problem has found applications 
in modern topics, such as jamming in granular media, colloids, and the structure of liquids 
and the glass transition [2]. Despite vast progress made in developing a statistical mechanics 
for such systems (3j, basic questions remain unanswered, including What is a 

jammed state and how to characterize its state of randomness? 

In an attempt to define the jammed states in a rigorous way, Torquato and coworkers 
have proposed three categories of jamming {(J: locally, collectively and strictly jammed. 
This problem is intimately related to the existence of well-defined upper and lower limits in 
the density of disordered packings; random close packing (RCP) and random loose packing 
(RLP) jj; a longstanding open question in the field. A definition of RCP requires proper 
definitions of jammed states and the concept of randomness These previous definitions 
of jamming are based purely on geometrical considerations, sufficient to describe frictionless 
grains 6] but not for granular materials where friction dominates Q. Figure UK illustrates 
the point: a frictionless hard sphere system is not locally jammed if only normal forces are 
considered, since the ball can freely move in the vertical direction. The same geometrical 
configuration is locally jammed if friction is allowed between the particles, revealing the 
importance of forces in the definition of jamming for frictional particles. 

The goal of the present paper is to characterize jamming and the degree of randomness 
for frictional hard spheres and obtain equations of state. Our framework is rooted in statis- 
tical mechanics and considers calculation of the entropy at the volume-force (V-F) ensemble 
level Js!] through force and torque balance conditions. The simulations and theory devel- 
oped allow understanding of the qualitative behavior of equations of state relating volume 
fraction, entropy and compactivity, and illuminate the nature of RCP and RLP. We follow 
three complementary approaches to characterize the entropy of jammed matter, (a) Com- 
puter studies: We first investigate frictional packings of spheres at the jamming transition 
numerically. We compute the equations of state, entropy and compactivity, as a function 
of volume fraction, ranging from RLP to RCP. The entropy is calculated by two numerical 
methods, direct analysis of volume fluctuations via a fluctuation-dissipation theorem and 
graph theoretical methods. Simulations show that random loose packings are more disor- 
dered and have higher compactivity than random close packings. In order to rationalize 
the simulation results we follow two theoretical approaches: (b) We develop the concept of 
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randomness in the V-F ensemble following the Gibbs distribution [3j , which is different from 
the measurement of randomness of a single packing in terms of an ensemble of order param- 
eters J5]. Since analytical expressions of the entropy are difficult to obtain in this ensemble 
we develop, (c), A mesoscopic ensemble approximation based on mechanical equilibrium 

than the minimum isostatic 

coordination as conjectured by Alexander |9j (see also 
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done under the mesoscopic approximations of 



12l . [13j ) . Calculations are 



giving rise to a mesoscopic configurational 



entropy, achieving a minimal value at RCP and maximal value at RLP. The results charac- 
terize the disorder of RCP and RLP at the mesoscopic level in general agreement with the 
simulations. They also suggest that the configurational entropy requires augmentation to 
include the entropy of the microscopic states neglected at the mesoscopic level. 

Numerical simulations. — We investigate computer generated packings of 10,000 spherical 
equal-size particles of 100/zm diameter interacting via Hertz (normal) and Mindlin (tangen- 
tial) contact forces (with shear modulus 29GPa and Poisson's ratio 0.2) with Coulomb 

nn 

friction, /i, using methods previously developed in [8L II 1|. Packings characterized by differ- 
ent fi are generated by compressing a gas of particles from an initial (unjammed) density, 
(pi, with a compression rate, T, until a final density, <p, at the jamming transition. 

The friction /i ranges from to 00 producing packings with coordination number varying 
from Z « 6 to Z « 4, respectively. Thus, the preparation protocol produces packings with 
densities parameterized by T, fa). In general, the lower F the smaller the obtained </>, 
while the larger fa the larger of the packing. We also find that under our numerical protocol 
there exits a common function Z(fi) over the different T and fa (see js]). For \i — > 00, <fi 
ranges from the RLP limit 0rlp ~ 0.55 obtained when r — > and 0j < 0.55 to the RCP 
limit 0rcp ~ 0.64 obtained for larger T and fa — ► 0.64 (note that merely changing T is not 
sufficient to allow 0rcp to approach 0.64, as discussed in [8J, for \i — > 00). For \i = 0, the 
density is approximately (f> ~ 0rcp for any (T,fa). For intermediate /x, the packings follow 
the phase diagram as obtained in Q]. 

The calculation of the entropy as a function of volume fraction is realized by using volume 



fluctuation analysis through a Fluctuation- Dissipation relation [14| complemented with 
graph theoretical methods Q|- We first define the Voronoi cell associated with each 
particle i and calculate its Voronoi volume Wj. We then perform statistical analysis of the 



volume fluctuations by considering a cluster of n contacting particles with volume W n 
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Y^Wi- We calculate the average volume, (W n ) and fluctuations ((W n — (W n )) 2 ), where (•) 
is an average over many n-clusters. We find that for sufficiently large n ~ 1000, contrasting 
the results of [17|, the fluctuations scale with n and therefore are extensive and well-defined 
(see inset of Fig. [Ib). 

From the large n behavior we extract the fluctuations, plotted in Fig [Tb, versus for 
every packing studied. The compactivity of the packing, X, is obtained via the integration of 
the fluctuation relation ((W„ - (W n )) 2 ) = XX 2 d(W n )/dX as X~ x = A d(W n )/((W n - 
(VV n )) 2 ), where we use that <fi(X — > oo) — > 0rlp [8j, and A is the analogue of the Boltzmann 
constant. Since Voronoi volumes are additive, (W n ) = (W) = NV g /(p, where V g is the 
volume of the grain. Therefore, the above integration is rewritten as: 

r 4>{x) 

(X/Vs)' 1 = A / #/(0 2 ((W„, - (W„» 2 », (1) 

where ((W n — (W n )) 2 is the fluctuation density. We may then utilize the fluctuation density 
as a function of 0, shown in Fig. [lb, and integrate along a line of constant Z(fi). We note 
that while the fluctuation densities for all Z(fi) in this study collapse onto a single curve, 
illustrated in Fig [lb, the limit of integration, 0rlp in Eq. (CD), changes as discussed in the 
phase diagram of [8|, increasing as /i decreases. 

The equation of state, <f)(X), is plotted in Fig. [lb for different values of the average 
coordination number of the packings, Z(fi), revealing that as we approach 0rcp ~ 0.64, 
X — > 0, regardless of the value of /i. Further, X — » oo as we approach 0rlp, with the 
smallest volume fraction of the RLP appearing for fi — > oo and Z w 4, with 0rlp ~ 0.55. 

The entropy, S, and its density, s = S/N, are obtained by integrating (X/Vg)' 1 = 
— 4> 2 ds/d(J), as: 

r4>nci> 

5(0) - s(0R CP ) = A / d0/[(X(0)/V 9 )0 2 ]. (2) 



This analysis provides the entropy up to a constant of integration s(0rcp) j . 



b obtain the 



entropy of RCP we use an independent method based on information theory [15j, [16j to obtain 
another estimation of the entropy. We use the Voronoi cell and Delaunay triangulation for 
each particle to define a Voronoi network. We construct a graph as a cluster of n contacting 
particles which, by means of graph automorphism, can be transformed into a standard form 
or "class" i of topologically equivalent graphs considered as a state with an occurrence p(i). 
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In practice, we determine p(i) by extracting a large number m of clusters of size n from the 
system and count the number of times, fc, a cluster i is observed, such that: p(i) = fi/m. 
The Shannon entropy of a clusters of size n is H(n) = — X^2p{i) lnp(z), and the entropy 
density is s = lim n _ +00 [i/(n + 1) — H(n)}, converging so rapidly that even moderate values 



of n are enough to obtain a sufficient approximation of s [16J . The Shannon entropy density 
provides an estimation of the entropy for the RCP state, s(0rcp) ~ 1.1A, serving as the 
constant of integration for the entropy density as realized by volume fluctuations. The 
resulting entropy density is plotted in Fig. [TH versus for different Z{p). 

When comparing all packings with different Z(p) and 0, the maximum entropy is at 
the minimum volume fraction of RLP 0rlp ~ 0.55 when X — > oo at infinite friction. The 
minimum entropy is found for the RCP state at </>rcp ~ 0.64 for X — > 0, for all the values of 
friction, indicating the degeneracy of the RCP state, contrary to the common belief that the 
RCP limit corresponds to a state with the highest number of configurations and therefore 
the highest entropy. 

Statistical mechanics of frictional hard spheres. — Next, we use the ensemble of jammed 
matter to rationalize the obtained equations of state 3J. Experiments of shaken grains, 
fluidized beds and oscillatory compression of grains [l4j indicate that granular materials 
show reversible behavior, and the analogue of the conserved energy, E, in thermal systems 
is the volume V = NV g /<p, for a system with iV at positions f£. Thus, the number of 
configurations, Q, and the entropy in the microcanonical ensemble of jammed hard spheres 
is defined as |3|] 

Q(V) = e s ^ x = J 8{V - W(r1))e jam (rl)2>l (3) 

Just as dE/dS = T is the temperature in equilibrium system, the "temperature" in 
granular matter is X = dV/dS. Here 0j am (^) is a constraint function restricting the 
integral to the ensemble of jammed states, W(r*j) is the volume function associated with 
each particle taking the role of the Hamiltonian. The crux of the matter is then to properly 
define Oj am and W to calculate the entropy and volume in the ensemble of jammed matter. 

Volume and force V-F- ensemble. — A minimum requirement of Ojam(^i) is to ensure touch- 
ing grains, and obedience to Newton's force and torque laws. As in the numerical simulations, 
the volume function, W(r*j), is taken as the volume of the Voronoi cell associated with each 
particle at position fj, for which an analytical form has been obtained in [8]. Following Eq. 
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FIG. 1: Numerical and theoretical results, (a) A disk in 2d under mechanical equilibrium by two 
nearest neighbor contacts. The disk is not jammed under a normal force interaction, but it jams 
when tangential forces are present, (b) (AW^) n versus (ft. Inset shows the extensivity of (AW^) 
versus n. (c) (ft versus X from the integration of (b). (d) Entropy versus (ft from the integration of 
(c). Entropy at RCP achieves a value of 1.1 A as calculated by the Shannon entropy at RCP. (e) 
Prediction of the mesoscopic theory for s me so(<^)/A + 1.1. (f) Prediction of the mesoscopic theory 
for (ft{X). 



PJ, the entropy in the V-F-ensemble of frictional hard spheres takes the form: 



VP 



D S{V)/\ 



S(y-W(fi))x 



II { 6 ( XJ fa) 6 ( ZJ fa x f v) 6 (fa - fa) x 



n[ (Mf-4 T )^([(^) 2 -i])% 



Vr, 



(4) 



where f. 



'j 



fj, the normal inter-particle force is f i 



N 



If, 



ij ' r ij\-, 



the tangential force: 



= \fij — (fa ■ fij)fij\. All quantities are assumed properly a-dimensional for simplicity of 
notation. The terms inside the brackets {■} correspond to the jamming constraint function 
Gj am in Eq. ([3]), and therefore define the ensemble of jammed states. The first three 
5— functions inside the big brackets impose Newton's second and third law. The Heaviside 
G— function imposes the Coulomb condition and the last 5— function imposes the touching 
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grain condition for hard spheres, assuming identical grains of unit diameter. Integration is 
over all forces and positions which are assumed to be equally probable as in the flat average 
assumption of the micro-canonical ensemble. 

The conditions specified in Eq. (J3J) are met in the numerical packings, thus the results 
of Figs. [Tb-d can be interpreted as the ensemble average of Eq. (j3J) under assumption of 
uniformity in the jammed configurations. However, Eq. (j3J) is difficult to solve. Analytical 
progress can be done by considering a coarse-graining of the Voronoi volume function and 
working with mesoscopic theory to obtain a configurational entropy at the mesoscopic 
level. 

Coordination number Z-ensemble. — Simple counting arguments, neglecting correlations 
between nearest neighbors, consider that a necessary condition for mechanical equilibrium 
is that the number of independent force variables must be larger or equal than the number 
of linear independent force/torque balance equations. Alexander Q conjectured that at 
the transition point for frictionless spherical packings 0, [l^] the system is exactly isostatic 
with a minimal coordination, Z = 2d = 6 in 3d. Such a conjecture can be extended to 

n 

the infinite friction case, where Z = d + 1 = 4 [13[ . In the presence of finite inter-particle 
friction coefficient /i, the analytic form of Z(/a) is difficult to achieve since the counting 
argument involves nonlinear inequality constraints through the Coulomb condition. Despite 
;he theoretical difficulty, there exists a dependency of Z with \i suggested by simulations 



12|. 



We consider a mesoscopic free volume function coarse-grained over a few coordination 

n 

shells that reduces the degrees of freedom to the coordination number z as shown in |8[: 
w(z') = ^jrVg. The w(z') function is calculated by obtaining the probability distribution 
of Voronoi volumes associated with a particle in the jammed state. This probability is 
decomposed into a bulk term depending on w and a contact term depending on z'. Note 
that z' refers to geometrical contacts that may carry no forces and can be larger than the 
mechanical coordination Z given by the isostatic condition: Z < z' jsj]. 

The mesoscopic approximation refers to the use of such a mesoscopic free volume function 
in the partition function, taking into account the effect of the environment of a particle rather 
than the particle itself. Thus, we can reduce the partition function to a single non-interacting 
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particle. The mesoscopic entropy density is obtained in the canonical ensemble [3j, 

s m eso = (w)/XX + A In j ^ g(z') exp (-^^) dz', (5) 

where the isostatic condition enters through the limit of integration. The space of config- 
urations is considered discrete since the states are collectively jammed p. If the typical 
distance between states is h z , the density of states is g(z') oc h% analogous to discretization 
of phase space imposed by the Heisenberg principle with density of states h~ d , where h is 
the Planck constant. Equation. (jSJ) is a coarse-grained mesoscopic form of the full entropy 
in Eq. (J4J) reducing the degrees of freedom to z. Thus the effective "Plank constant" h z 
appears in Eq. (jSJ) and not in Eq. (jlj). 

The mesoscopic entropy of Eq. (jSJ) is plotted in Fig. [Tfe as a function of 4> for different 
values of Z. We see that it captures the general behavior found in the simulations, i.e., 
maximal at RLP for Z = 4 and minimal at RCP. Furthermore, all the curves for different 
Z approach S ~ lnX as X — > 0, similar to a thermal ideal gas. We conclude that, at 
the mesoscopic level, the entropy vanishes at RCP (in fact it diverges to — oo when <fi — > 
0rcp closer than a constant proportional to h z ), providing a characterization of RCP. This 
result qualitatively resembles behavior of the complexity of the jammed state in the replica 
approach to jamming [18| . We use h z = 0.01 in Fig[Tfe such that the mesoscopic entropy 
vanishes very close to the predicted value of 0rcp ~ 0.634 [8| and the value of the microscopic 
entropy as calculated via graph theoretical methods. The value of h z is chosen to fit the 
theory with simulation as close as possible, where the only constraint imposed by theory is 
h z < 1. An important result is the direct implication of a larger number of states available 
to jammed systems at RLP with respect to any higher volume fraction. 

It is possible to interpret the RCP as a Kauzmann point (K-point) [3| in analogy with 
the density (temperature) at which the configurational entropy of a colloidal (molecular) 
glass vanishes at the ideal glass transition. 

We augment the mesoscopic entropy with the entropy of the microscopic states to obtain 
the full entropy as s = s mcso + s m i cro . Since for frictionless packings ,s m eso vanishes, then 
s = s m i cr o, implying that we can obtain s m i cro from the full entropy of the K-point as 
calculated numerically. Therefore, Fig. [Tfe is plotted as s meso /X + 1.1 to obtain a plot of 
the full entropy to compare with the full ensemble Eq. (J3J) and simulations. This constant 
is obtained from the value given by the Shannon Entropy, regarded as the entropy of RCP, 
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which is 1.1A. 

In summary, a notion of jamming is presented that applies to frictional hard spheres, as 
well as frictionless ones. The entropy reveals interesting features of the RCP and RLP states 
such as the fact that RLP is maximally random with respect to RCP and that both limits 
can be seen in terms of the entropy and equation of state. The theoretical model captures 
the shape of the entropy and the equation of state, but not the actual values at RLP. Future 
study will be devoted to capture the volume fluctuations not only at the mesoscopic level 
but also at the microscopic level neglected in the present theory, which may lead to an 
understanding of the finite value of entropy found at RCP in the numerics. The agreement 
between theory and simulation is sufficient to indicate that the methods presented herein 
are appropriate for evaluating the entropy of jammed matter allowing the characterization 
of the state of randomness of RLP and RCP. We note that an interesting recent work [3] 
finds random very loose packings with negative effective temperature in an energy ensemble 
approach. Such states (analogous to negative X in the volume ensemble) are discussed 
elsewhere js|. Our results are in general agreement with those of {^(J. 
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